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ABSTRACT 

We consider a situation in which a pulsar is formed inside or close to a high density region 
of a molecular cloud. Right after birth, the pulsar was very active and accelerated hadrons 
and leptons to very high energies. Hadrons diffuse through the supernova remnant (SNR) and 
some of them are trapped in the nearby cloud interacting with the matter We extend a recent 
time-dependent model for the 7-radiation of pulsar wind nebulae (PWNe) to describe this 
more complicated astrophysical scenario. The example calculations have been performed for 
two objects, IC443 and W41, which have recently been discovered as sources of TeV 7-rays. 
In this model the low energy TeV emission should be correlated with the birth place of the 
pulsar and the region of dense soft radiation rather than with its present position, provided that 
the injection rate of relativistic particles into the PWNa has been much more efficient at early 
times. The high energy TeV emission should be correlated with the location of dense clouds 
which were able to capture high energy hadrons due to their strong magnetic fields. 

Key words: gamma rays: theory — supernovae remnants — ISM;individual (IC 443, W41) 
— pulsars 



1 INTRODUCTION 

A few recently discovered TeV 7-ray sources seem to be clearly 
related to pulsar wind nebulae (PWNe) 1 Aharonian et al. 2006^). 
However, in some cases the TeV source is clearly displaced from 
the present position of the pulsar. This can be explained either by 
the fact that while the pulsar was very young it injected particles 
and then it moved significantly during its lifetime of ~ 10* years or 
by the displacement of the low energy pairs from their place of ori- 
gin due to the interaction of the PWNa with the reverse shock of the 
supernov a remnant (SNR) or b y a back flow in the bow shock neb- 
ulae (e.g. iBlondin et all feOOlh ). All these explanations are likely 
and they cannot be excluded at the present stage of knowledge. 
Here, we consider the first scenario in more detail applying it 
to two such example TeV 7-ray sources with eviden ces of asso- 
ciated energetic pulsars: IC 443 (MAGIC J06 16+225, [Albert et all 
1 I2OO7I), see also Humensky et al. (2007 )) and W41 ( HESS J1834- 
OSV- lAharoniarTet al.l ( l2006ah and I Albert et alj ( l2006b ). 

In the case of IC 443 the VHE 7-ray source is located in 
the direction of a dense molecular cloud which probably lies in 
front of the parent asymmetric SNR. Moreover, an X-ray nebula, 
G 189.22+2.90, has been discovered near the s hell of IC443 by re- 
cent XMM-Newton and C handra observations jBocchino & Bvkovl 
I200ll : l01bertetal.ll200lh . It is interpreted as being due to the pres- 
ence of an energetic Vela-type pulsar. The X-ray PWNa is displace d 
from the location of the TeV source by ~ 20', jAlbert et al.ll2007h . 
We suggest that this pulsar can also be responsible for the observed 
TeV 7-ray source. The 7-rays can be produced by particles accel- 
erated inside the SNR during the early age of the pulsar. In fact, a 



pulsar with a velocity of ~250 km/s should change its position by 
^ 18' (at a distance of 1.5 kpc) by the time it reaches the age of 
3 X 10'' years. 

The VHE 7- ray source HESS J1834-087 looks extende d with 
a size of ~ 12' jAlbert et al.]|2006l : lAharonian et alj|2006ah . It is 
positionally coincident with the shell-type SNR, G23.3-0.3 (W41), 
with an age estimate of ~ 8 x 10"* years. A Vela-type pulsar, PSR 
J 1833-827, with the period of 85 ms, is also p ossibly connected 
with this supernova l iGaensler&Johnstoiilll995h . Its characteristic 
age and distance are consistent with the age and distance of the 
SNR W41. However, the pulsar is displaced from the center of the 
SNR (and also from the position of the TeV source) by ~ 24'. Such 
positional disagreement might again be explained by the movement 
of the pulsar during ~ 10^ years with a velocity of ~ 250 km/s. 
The pulsar seems to be located in the direction of one of the exten- 
sions of the TeV 7-ray source ( lAharonian et al ]|2006al) . 

In this paper we extend the hybrid (leptonic-hadronic) time 
dependent radiation model for PWNe (Bednarek & Bartosik 2003, 
hereafter BB 03) for the more complicated case of fast-moving pul- 
sars in the vicinity of dense molecular clouds. 



2 A MODEL FOR 7-RAY PRODUCTION 

We consider a scenario in which an energetic, fast-moving pulsar is 
created in a supernova explosion close to a dense cloud. The evolu- 
tion of the Pulsar Wind Nebula, created by such a pulsar, can differ 
significantly from the case of an isolated PWNa of the Crab type. 
The higher density of the surrounding medium results in a slower 
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expansion of the SNR and lower adiabatic losses of particles. The 
high density matter can provide additional soft radiation and mat- 
ter targets for the particles accelerated in the vicinity of the pulsar. 
Moreover, the pulsar, changing position in time, injects relativistic 
particles (lept ons and hadrons) in different places. 

The IBBI losh radiation model is a time-dependent, but one- 
space-dimensional model for isolated PWNe. In order to consider 
in detai l the time-dependent radiation processes inside the neb- 
ula |bb| losh consider a rather simple model for the expansion of 
the PWNa and th e SNR based on the early models proposed in 
the literatu re (e.g. tostriker & Gunn ' 197 iMPacini & Salvati 1 19731 : 
iRees & G unn 1974; Reynolds & Chevali erTr984h y Note however, 
that it is at present well known that the evolu tion of such systems 
is quite complicated (as considered by e.g. liBlondin et al.|[200ll : 



Iv. d. Swaluw et al.ll200ll : iBucciantini et al. l2003h ). At the present 
stage we are not able to take into account the complicated evolu- 
tion of the nebula envisaged in these papers. It can show significant 
differences in the case of nebulae with an age larger than a few 
thousand years. Therefore, we limit our considerations to a sim- 
ple evolution model of the nebula and concentrate on a detailed 
description of the complicated time-dependent radiation processes 
inside the nebula. As we show farther, our model considers the ac- 
celeration of leptons at the relatively earlier stage of the PWN evo- 
lution. Therefore, the exact evolution of the nebula at later stages 
does not have a very important effect on the radiation processes. 
Still, it might have an effect on the morphology of the 7-ray source. 

Ou r model applies the h y pothesis of Arons and c ollaborators 
(see e.g. lHoshino et alj ( ll992h : lGallant & Arorii ( Il994b ) according 
to which leptons gain energy by being resonantly scattered off the 
Alfven waves generated above the pulsar wind shock by hadrons 
injected by the pulsar. The model assumes that most of the rota- 
tional energy loss rate of the pulsar is taken by relativistic hadrons 
which can gain energy corresponding to approximately 20% of 
the maximum potential drop through the pulsar magnetosphere. 
Therefore , hadrons are injected mono-energetically. According to 
iHoshino et al.l l[l992i) leptons gain a power-law spectrum with an 
index not far from two. In our model we assume this power-law 
index to be independent of time. 

These leptons and hadrons diffuse inside the SNR and loose 
energy in adiabatic and radiation processes. The diffusion of 
leptons is much slower than the one of the hadrons due to 
huge radiation losses on synchrotron, bremsstrahlung and IC 
processes. Therefore, hadrons are able to partially escape from 
the SNR, contrary to leptons. Note, however, that we apply a 
rath er simple model for the diff usion of hadrons from the neb- 
ula teednarek & Protheroell2002h . This model does not take into 
account the Rayleigh-Taylor instabilities on the border separating 
the PWNa from the SNR. Due to these instabilities, the escape of 
hadrons from the nebula may become more efficient and the tar- 
get for relativistic hadrons which are still inside the PWNa may 
increase. These two effects have an opposite influence on the 7-ray 
spectrum produced by hadrons inside and outside the PWNa and 
are not taken into account in our model. 

Due to the evolution of the SNR in time (radius, expansion 
v eloci t y , density of matter, magnetic field strength, see Eqs. (1-7) 
in lBBi losl) ) the above mentioned processes are time-dependent in a 
non-trivial way. The evolution of the equilibrium spectra of leptons 
is described by E qs. (11-12) of B B (03 ) and the one of hadrons 
by Eqs. (10-11) of iBednarek & BartosikI j2004l) . In the considered 
time-dependent scenario these equations cannot be solved analyt- 
ically. We apply a n iterative procedure, the time-step method, de- 
scribed inHgl lol). 



Based on the equilibrium spectra of particles we calculate 
the 7-ray production from the SNR in different radiation pro- 
cesses (bremsstrahlung and inverse Compton (IC) for leptons and 
interactions with matter for hadrons). Hadrons which escaped 
from the SNR can be partially captured by molecular clouds and 
also produce 7-rays. In the calculations of the 7-ray spectra for 
hadro ns we apply the scaling break m odel for hadronic interac- 
tions dWdowczvk & Wolfendafll 19871) suitable for the considered 
energy range of hadrons. 

In order to be able to apply this time-dependent model for 7- 
ray production from PWN to the considered astrophysical situation, 
it has to be extended by two issues: the fast motion of the pulsar and 
the vicinity of a dense molecular cloud. 

The analysis of the radiation processes around pulsars is 
considerably complicated in the case of pulsars which are able 
to significantly change their later position with respect to their 
birth place. This situation obviously concerns the Vela-type pulsars 
which reach large velocities at birth and which are at later evolution 
stages still energetic enough to produce magneto- spheric 7-rays (7- 
ray pulsars). In their case, the radiation processes described above 
can be additionally distributed along the path of the moving pulsar. 
In such a complicated case, the 7-ray source can be shifted with re- 
spect to the present position of the pulsar depending on the history 
of the injection of relativistic particles and the distribution of the 
target (soft radiation, matter) close to the path of the pulsar. 

In order to compute the spectrum and the morphology of 
the VHE 7-ray source a full three-space-dimensional and time- 
dependent treatment is necessary. Unfortunately, this is too diffi- 
cult (requires consideration of three-dimensional radiation and dif- 
fusion processes) to be constructed in detail at present. We limit 
our calculations to a spherically symmetric case, which can give a 
correct estimate for the 7-ray spectrum and flux since the soft radi- 
ation fields (CMB, infrared and optical background photons in the 
Galactic plane or produced by the molecular cloud) do not depend 
on the movement of the pulsar. Note, that the IC scattering of the 
synchrotron radiation is negligible for old PWN. However, in such 
an approach we are not able to calculate the 7-ray source morphol- 
ogy. In this spherically symmetric model, we can take into account 
the additional soft radiation target, created by the emission from 
the molecular cloud, but we cannot take into account the possible 
influence of the cloud on the spherical geometry of the expanding 
SNR and the related PWNa. Therefore, this model can give an ap- 
proximately correct description of the considered picture only in 
the case that the cloud is not too close to the expanding SNR. 

There are arguments (see Sect. 3 in IBBI l losh ) that the injec- 
tion rate of particles into the nebula can be limited to a specific 
period after the pulsar formation. The processes occurring close to 
the pulsar surface or inside its inner magnetosphere are only ac- 
tive i.e. during the first 10"^ — 1 0^ years after the pulsar formation, 
see e.g. lUsov & Melro"s3 ( Il995[) . For example, nuclei, which are re- 
sponsible for the acceleration of leptons, are only injected when the 
surface temperature of the neutron star is high enough. 

In our calculations we apply the basic parameters of the ex- 
ploding supernova as used in our earlier modeling of the Crab Neb- 
ula (BB 03), i.e the initial period of the pulsar is 10 ms, the initial 
velocity of the supernova shock is 2000 km s~^, and the mass of 
the nebula is 4 Mq. We assume that the supernova exploded in a 
medium with a density of ~ 20 part, cm^^ (Dame et al. 1986) 
and a magnetic field of 10~^ G. Such a large density is supported 
by the presence of dens e clouds around the SNR . Following Arons 
and collaborators (e.g. iGaUant & Aron3 ( Il994l) ). we assume that 



hadrons take 95% of the rotational energy loss of the pulsar. Only 
5% of the energy of hadrons is transferred to leptons. 

The infrared soft radiation field is expected to be considerably 
stronger inside the SNR due to the emission from the dust heated 
by the supernova shock. Following other works, we apply the soft 
radiation field of the interstellar medium at the galactic disk as de- 
fined in toe Jager & Mastichiadi3ll997l : iGaisser et al.lll998h . It is 
composed of the cosmic microwave background radiation (CMB) 
with a temperature of 2.7 K, the infrared background with an en- 
ergy density two times larger than the CMB and a temperature of 
25 K, the optical background defined by energy densities equal to 
the CMB and with characteristic temperatures between 5000 K and 
IC* K. Moreover, we add the infrared radiation field with an energy 
density three times larger than the CMB and with a temperature of 
45 K, which is supposed to be pro duced by dust inside the molecu- 
lar cloud (based on the IRAS data jMufsonetalll986h ). 
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Table 1. Parameters of the model considered for IC 443 and W41 





IC443 


W41 


pulsar or PWNa 


G189.22+2.90 


PSRJ1833-827 


pulsar surface magnetic field 


3 X 10^2 G 


1.2 X 1012 G 


present pulsar period 


145 ms 
(estimated) 


85 ms 
(observed) 


age of the pulsar 


3 X lO"* yrs 


8 X lO"* yrs 


duration of lepton injection phase 


10'' yrs 


4 X 10'' yrs 


spectral index of leptons 


2.4 


2.4 


density of suiTounding medium 


20 cm-» 


20 cm-3 



magnetic field of surrounding medium 10 ^ G 10 G 



3 IC 443/MAGIC J0616+225 

The parameters of the pulsar which creates the PWNa observed 
close to the SNR IC 443, G 18 9.22+2.90, are unk nown. But from 
comparison with other nebulae, [oibert et ^ ( 1200 ih argued that the 
pulsar has to be of a Vela-type with a present period of ~145 ms, 
a surface magn etic field of ^SxlO '^^ G, and an age of ~3 x 10^ 
years, see also ISocchino & Bvkovl /2001). The displacement be- 
tween the present location of the PWNa and the center of the 
SNR can be explained by the motion of the pulsar with a veloc- 
ity of ~250 km/s taking into account its estimated age. IC 443 
is an asymmetric shell-ty pe S NR with a diameter of 45' at a dis- 
tance of about 1.5 kpc ( lFesenlll984l: IClausen et al. 199 ^. It is a 

S iromi n ent X-ray source, see e.g. iTroia et alj ( i2006) : .Oibert et al.l 
200 ih : iGaensler et d] bOOd) . The EGRET has detected a 7-ray 
source above 100 MeV in the IC 443 SNR, 3EG J06 17+223 8 
(lEposito et al. 1996; Hartman et al. 1999). There is a large amount 
of molecular mass (< 10^ Mp) ) consistent with the distance to the 
SNR IC 443 tTorres et al.l2003h . The highest CO intensity detected 
is directly superimposed on the central positio n of the VHE 7-ray 
source MAGIC J06 16+225 (Albert et al '2007^. 

The pulsar with the parameters mentioned above resembles 
the 7-ray pulsars in Vela and PSR 1706-44. In fact, when we com- 
pare the 7-ray spectrum observed from e.g. PSR 1706-44 with the 
spectrum of the EGRET source 3EG J0617+2238, which is ob- 
served in the direction of IC 443, we get good consistency of the 
spectral shapes and fluxes (within a factor of two) of these two 
sources (see Fig. 1). Note, that the distances of the pulsars have 
significant errors due to non-exactly known parameters of the inter- 
stellar space. Therefore, we propose that, in spite of only marginal 
directional agreement, the EGRET source 3EG J0617+2238 is due 
to the pulsar responsible for the PWNa observed close to IC 443. 

To calculate the 7-ray spectra produced by leptons in the IC 
and bremsstrahlung processes we used the above described model 
with parameters of the supernova remnant, the pulsar, and sur- 
rounding medium mentioned at the end of section 2 and in Tab. 
1 . We investigated the dependence of the flux level and spectra on 
the duration of the lepton injection phase and their spectral index, 
keeping all other parameters fixed. We found that the MAGIC data 
can be fitted well by the values of these parameters reported in Tab. 
1. The difference between the present location of the PWNa and 
the position of the MAGIC TeV 7-ray source can be explained by 
assuming that leptons are accelerated only during the early phase 
of the PWNa. 



The population of high energy electrons, which produces the 
VHE 7-rays by IC scattering and bremsstrahlung, also produces 
synchrotron radiation in the magnetic field of the SNR. However, 
these old electrons do not have a sufficiently high enough en- 
ergy to produce synchrotron X-rays. The characteristic energy of 
synchrotron photons can be estimated by esynch ~ me7e-B/(4 • 
10^"^ G). Therefore, 10 TeV electrons would produce synchrotron 
emission with a characteristic energy of only about 50 eV. 

We calculate the 7-ray spectrum from the decay of 7r"s pro- 
duced by hadrons accelerated by the PWNa and interacting with 
the matter inside this nebula. However, since the nebula is rela- 
tively old, most of the hadrons accelerated during the early stage of 
the nebula have already escaped from it. Therefore, the level of the 
7-ray emission produced by hadrons which are still inside the SNR 
is relatively low (see Fig. 1). 

We also estimate the 7-ray fluxes produced by those hadrons, 
which escaped from the PWNa in the past activity stage of the pul- 
sar. It is very difficult to reliably estimate the part of hadrons which 
are captured inside dense clouds, since this process depends on the 
geometry of the cloud and the nebula, their relative distance, and 
above all on the geometry and distribution of the magnetic field in 
the region surrounding the PWNa. This accumulation of hadrons 
may (or may not) also depend on the energy of the accelerated 
hadrons, depending on the details of the capturing process (e.g. 
whether it is mainly due to energy-dependent diffusion or due to 
advection from the wind from the pulsar). Since we are not able 
to take into account these complicated processes, it is assumed 
that a small part of all hadrons escaping from the PWNa is cap- 
tured inside dense clouds. The example calculations are performed 
for the parameters of the PWNa related to IC 443, assuming the 
presence of a nearby cloud with a density of 10"^ particles cm^^ 
and an efficiency of accumulation of hadrons equal to 10^^. It 
is clear that, even for such a very low capturing efficiency, these 
hadrons can produce large fluxes of 7-rays in the energy range 
of ~1 0-10^ TeV. Note, that the upper limits to the 7-ray fluxes 
above 137/200 TeV by the Cygnus / CASA-MIA collaborations 
dAllen et al.ll 19951 : iBorione et al.lll99a) are more than one order of 
magnitude above our predicted flux. These upper limits constrain 
the capturing factor of hadrons in the molecular cloud to be less 
than 10"". 
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Figure 1. The 7-ray spectra (SED) produced by leptons in the IC (blue full curve) and bremsstrahlung (red dotted curve) processes inside the SNRs IC 443 
(left figure) and W41 (right figure). The 7-ray spectra produced by hadrons inside the supernova remnant (green dot-dashed curve), and by hadrons which 
escaped from the supernova remnant and which were captured by nearby dense clouds (black dashed curves). We also show the 7-ray spectrum observed from 
the inner magnetosphere of the Vela-type pulsar (PSR 1706-44) by the EGRET telescope (green triangles), unsealed for IC 443 due to comparable distances 
and scaled to take into account the different distances of PSR 1706-44 and W41 . The EGRET data on IC 443 are portrayed by black stars. The MAGIC VHE 
7-ray data are portrayed by blue points. The cyan and magenta an'ows indicate 90% upper limits to the UHE 7-ray flux from the CASA-MIA and CYGNUS 
collaborations jAUen et alJl995l ; lBorione et aljl995h . The parameters of the model are given in table[T] 



4 W41/PSR J1833-827/HESS J1834-087 

The pulsar, PSR J 1833-827, has been observed 24' away from the 
center of the TeV 7-ray source HESS J1834-087. The pulsar has a 
period of 85 ms (belongs to the Vela-type pulsar population), has 
an estimated sur face magnetic field of ~ 1 0^^ G, and lies at a dis- 
tance of ~4 kpc dOaensler&JohnstorJ I I995I) . The characteristic age 
of the pulsar, ~ 10^ years, and distance are consistent with the cor- 
responding par ameters of the S NR G23. 3-0.3 (W41, i.e. an age of 
8 X 10 years i lTian et aLllioOTl ). The observed displacement be- 
tween the present position of the PSR J1833-827 and the center of 
W41 can be easily explained by applying the above-mentioned age 
of the pulsar (and SNR) and assu ming that the v elocity of the pulsar 
is equal to 250 km s^^. Recentlv, iLandi et alj (12006) found a faint 
X-ray source within th e area of HESS J1 834-087AV41 in data from 
the Swift satellite and iTian et al. I j2007h found an extended X-ray 
feature spatially coincident with the VHE 7-ray emission. As in the 
case of IC 443, W41 is associated with a large molecular complex 
called "[23,78]" jPame et aljigs^) at a similar distance. 

We also apply the time-dependent model for the acceleration 
of hadrons and leptons in the vicinity of energetic pulsars for this 
object. The general parameters of the medium, the target photon 
fields with which particles interact, the initial parameters of the su- 
pernova remnant and the pulsar are kept the same, as in the case 
of IC 443. We found that the MAGIC data can be fitted well by 
the same spectral index of leptons as in the case of IC 443, but 
by a four times larger duration of the lepton injection phase. The 
results of calculations are also shown in Fig. 1. In our interpreta- 
tion, the differences between the 7-ray spectra produced by these 
two objects are mainly due to different values of the pulsar sur- 
face magnetic fields and the duration of the lepton injection phase. 
They determine the energy distribution and the power of the ac- 
celerated particles as well as their energy losses in case of leptons 
(synchrotron and indirect IC process) and the diffusion of hadrons 
inside the nebula. Note, that the leptons produce 7-rays not at the 
present position of the pulsar, but closer to its birth position and the 
position of the molecular cloud due to the stronger infrared photon 
field of the cloud. 

We also calculate the 7-ray spectra from the decay of pi- 
ons produced by hadrons, accelerated in the vicinity of the pulsar. 



which are still inside the SNR, and those ones which escaped from 
the SNR and interact with dense clouds. 

Note, that we scaled the EGRET data of PSR 1706-44 in the 
right plot of Fig. 1 to account for the different distances of PSR 
1706-44 and W41. In this way we predict the expected level of the 
flux in the EGRET energy range from the pulsar inside W41 which 
will be probed by the AGILE and GLAST telescopes. 



5 CONCLUSION 

We interpret the recently detected VHE 7-ray emission from two 
objects, IC 443 and W41, in terms of the time-dependent leptonic- 
hadronic model for the radiation proces ses i n SNRs containing 
Vela-type pulsars. We modify the model (BB 03) in order to take 
into account the presence of a nearby high density region which can 
serve as a target for relativistic hadrons, escaping from the SNR, 
and as an additional source of infrared photons upscattered to 7- 
ray energies by leptons inside the SNR. We show that the higher 
density of the surrounding medium results in a slower expansion 
of the SNR and lower adiabatic losses of particles in respect to 
isolated SNRs. Moreover, the nearby high-density matter acts as 
a target for hadronic interactions and provides additional infrared 
radiation fields for more efficient leptonic interactions. All of this 
causes higher 7-ray fluxes. 

We show that, in the case of IC 443 and W41, the observed 7- 
radiation can be explained by such SNR-PWN-pulsar systems for 
likely initial parameters of the pulsars and the SNRs and the sur- 
rounding medium. We suggest that the pulsar, the 7-ray sources at 
different energies, and the center of the parent SNR do not need 
to be co-spatial in the case of fast-moving pulsars close to a high- 
density region when the acceleration of leptons is much more effi- 
cient at the early age of the pulsar. 7-ray sources at energies below 
a few TeV are likely to appear in the place where the pulsars have 
been bom and the 7-ray sources at higher energies should be related 
to the location of the massive nearby clouds. On the other hand, the 
GeV 7-ray source should be co-spatial to the present position of the 
pulsar. 



Gamma-ray emission from PWNe 5 

ACKNOWLEDGEMENTS 

This work is supported by the German MPG and the Pohsh MNil 
grant 1P03D01028. 



REFERENCES 

Aharonian, F. A. et al., 2006a, ApJ 636, 777 
Albert, J. et al., 2006, ApJ 643, L53 
Albert, J. et al., 2007, ApJ, 664, L87. 
Allen, G. E. et al., 1995, ApJ, 448, L25. 

Bednarek, W., Protheroe, R. J., 2002, Astropart. Phys., 16, 397. 
Bednarek, W., Bartosik, M., 2003, A&A 405, 689 (BB03). 
Bednarek, W., Bartosik, M., 2004, A&A 423, 405 
Blondin, J. M., Chavalier, R. A., Frierson, D. M., 2001, ApJ, 563, 
806. 

Bocchino, F, Bykov, A., 2001, A&A 376, 248 

Borione, A. et al., 1995, Proc. 24th ICRC (Rome), 2, 439. 

Bucciantini, N., Blondin, J. M., Del Zanna, L., Amato, E., 2003, 

A&A 405, 617. 
Clausen, M. J., et al., 1997 ApJ 489, 143. 
Dame, T. M. et al., 1986, ApJ 305, 892 
De Jager, O. C, Mastichiadis, A., 1997, ApJ 482, 874 
Esposito, J. A. et al., 1996, ApJ, 461, 820. 
Fesen, R. A., 1984, ApJ 281, 658 

Gaensler, B. M. & Johnston, S., 1995, MNRAS, 275, L73. 

Gaensler, B. M. et al., 2006, ApJ, 648, 1037. 

Gaisser, T. K., Protheore, R. J., Stanev, T, 1998, ApJ 492, 219 

Gallant, Y. A., Arons, J., 1994, ApJ 435, 230 

Hartman, R. C. et al., 1999, ApJS 123, 79 

Hoshino, M. et al., 1992, ApJ 390, 454 

Humensky, T. B. (VERITAS collaboration), 2007, Proc. 30th 

ICRC, Merida, arXiv:0709.4298 
Landi, R. et al., 2006, ApJ, 651, 190. 
Mufson, S. L. et al., 1986, AJ 92, 1349 
Olbert, C. M. et al., 2001, ApJ 554, L205 
Ostriker, J.R, Gunn, J.E. 1971, ApJ 164, L95 
Pacini, F, Salvati, M. 1973, ApJ 186, 249 
Rees, M.J., Gunn, J.E. 1974, MNRAS 167, 1 
Reynolds, S.R, Chevalier, R.A. 1984, ApJ 278, 630 
van der Swaluw, E., Achterberg, A., Gallant, Y.A., Toth, G., 2001, 

A&A, 380, 309. 
Tian, W. W. et al., 2007, ApJ 657, L25 
Torres, D. F et al., 2003, Phys. Rept., 382, 303. 
Troja, E., Bocchino, F & Reale, F, 2006, ApJ, 649, 258. 
Usov, V.V., Melrose, D.B., 1995, Austr. J. Phys. 48, 571 
Wdowczyk, J., Wolfendale, A. W., 1987, J.Phys. G 13, 411 



